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Abstract 

We study two models describing the interaction of a two-level system with two quantum field 
modes. The first one is equivalent to a dissipative two-state system with just two boson fields in the 
absence of tunneling. The second describes two orthogonal fields interacting with the corresponding 
orthogonal dipoles of a two-level system. We show that both models present a partial two-mode 
SU{2) symmetry and that they can be solved in the exceptional case of resonant fields. We study 
their ground state configurations, that is, we find the quantum precursors of the corresponding 
semi-classical phase transitions, as well as their whole spectra to infer their integrability. We show 
that the first model in the exceptional case is isomorphic with the quantum Rabi model and allows 
just two ground state conhgurations, vacuum and non-vacuum. The second model allows four 
ground state configurations, one vacuum, two non-vacuum single mode and one non-vacuum dual 
mode, and give analytic and numerical pointers that may suggest its integrability. We also show 
that in the single excitation subspace these models can serve as a fast SU (2) beam splitter even in 
the ultra-strong coupling regime. 
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I. INTRODUCTION 


The so-called quantum Rabi model [T], 

Hr = + uja^a + g{d^ + d) (1) 

modeling the interaction of a two-level system, described by the transition frequency ujq and 
the Pauli operators dj with j = x, y, z, with a boson held, described by the held frequency 
u and the annihilation (creation) operators d (al) can be seen as a single-held version of the 
dissipative two-state system [2] in the absence of tunneling, A = 0, 

Hr = - yda, ^ (^jd]dj + ^ Qj (^d] + dj^ d^. (2) 

j j 

The dissipative two-state model is characterized by a spectral function, J{oj) = vr gj6{u}j — 
u) and is solvable, for example, for sub-Ohmic, Ohmic and super-Ohmic spectral functions, 
J{uj) oc ca® with s < 1, s = 1 and s > 1, in that order [2]. On the other hand, the solvability 
and integrability of the quantum Rabi model has been recently discussed for any given 
parameter set [3H8]. An equivalent approach has been used to explore the integrability and 
exceptional solutions of the two-qubit quantum Rabi model [nHS]. 

Here, we are interested in exploring the exceptional solutions of models describing a 
single qubit coupled to just two boson helds. Our motivation is twofold. First, the use of 
Bargmann approach may render the system integrable for exceptional parameter sets in a 
way equivalent to that found in the two-qubit single-held case. Second, circuit quantum 
electrodynamics (circuit-QED) may provide a direct testing ground for such a model from 
weak to ultrastrong couplings [IT] . 

2 2 

Hi = Y^z + Ujd]dj + Y 9j (aj + (3) 

i=i i=i 

and cavity-QED may provide an equivalent model by Raman adiabatic driving of a four-level 
atom coupled to two cavity electromagnetic held modes [15] . 

2 

H2 = Y^^ + + 9l («1 + Olj da: + ig2 (dl - 02 j dy. (4) 

i=i 

Note that both models conserve parity. If = dehned in terms of the total number of 

excitation N = d^j^ + d\di + dld 2 + 1/2. 
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This manuscript is structured as follows. First, we will study possible candidates for 
symmetries on these models and regimes where they are equivalent to well known models. 
Here, we will introduce the exceptional case provided by resonant helds, where Hamiltonian 
Hi reduces to the standard quantum Rabi model and Hamiltonian H 2 is invariant to a 
SU{2)^SU (2) transformation for identical couplings. Then, we will focus on this exceptional 
case of resonant fields to study the ground state structure of both models. While Hamiltonian 
Hi shows a simple ground state conhguration that includes just a vacuum and a non¬ 
vacuum ground states, Hamiltonian H 2 shows a more interesting ground state conhguration 
landscape with four possible conhgurations, one vacuum, two non-vacuum single mode, and 
one non-vacuum dual mode. Next, we will discuss the integrability of Hamiltonian Hi due 
to the isomorphism with the quantum Rabi model and give analytic and numeric arguments 
that point in the same direction for Hamiltonian H 2 . Finally, we will demonstrate that the 
partial SU{2) symmetry, shown by both models in the exceptional case of resonant helds, 
allows us to construct closed form evolution operators in the weak coupling regime. We 
will use these evolution operators to show that these models may be used as SU{2) beam 
splitters in the single excitation subspace and, hnally, we will present a couple of exact 
numeric time evolution that show this beam splitter ehect survives even in the ultra-strong 
coupling regime for short evolution times. 

II. SYMMETRIES AND EQUIVALENCE WITH OTHER MODELS. 

First, we want to bring forward that both Hamiltonians Hi and H 2 are invariant to full 
rotations, 6 = 2mi with n = 0 , 1 , 2 ,..., under the unitary transformation, 

in other words, 

U(2u7r)i7jf/^(2u7r) = H^. (6) 

The held part of this transformation is related to Schwinger two-mode representation of 
SU( 2 ) [in], J+ = 0 ^ 02 , J- = 0102 , Jo = (oiOi — 0202 ) 72 . This does not provide us with any 
information but note that in the case of identical qubit-held couplings, gi = g 2 , Hamiltonian 
H 2 is invariant to any given rotation parameter [T5] . 

U{e)H2U\e) = H 2 , gi = g2. (7) 
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On the other hand, it is well known that there exists an exact nnitary transformation that 
maps the dissipative two-level model, Hl, into a linear nearest neighbor chain of coupled 
bosonic modes where just the hrst one of them is coupled to the qubit mi. It is not surprising 
that such transformation in the hnite case is related to the SU(2) unitary displacement 
operator, 

D(0 = tane = -• (8) 

9i 

This two-mode displacement yields an effective model where the qubit couples to only the 
hrst boson held in the usual quantum Rabi model form, and the hrst and second boson held 
couple between them with a beam splitter form, 

H,n = DiOHiD^iO, ( 9 ) 

2 

= ^ ^ja]aj + X (^a\a2 + dia^ j -t- S' Si j (10) 

i=i 

Here, we have dehned effective held frequencies fli = (cuiSi + '^25'!) 19^ + '^ 2 fi'i) /s'^, 

effective held coupling constant A = { 0 J 2 — oJi) 9192 /9“^, and effective qubit-held coupling 
9 = + 92 - Note that choosing tan,^ = — 5 'i/fi '2 as two-mode displacement parameter 

just interchanges the boson held modes. 

A set of unitary transformations cannot bring Hamiltonian H 2 into an expression similar 
to Hid, thus we are reduced to explore regimes where they may be equivalent. We can start 
by using the same two-mode displacement on H 2 and hnd. 


H 2 d = m)H 2 D\i), 


( 11 ) 


2 


d'z + E ^jdijOij -f- A ^s|ci2 T -|- 


i=i 

, 9l- 9l^ 2^W2. 

SSi -|--Si —-S2 

9 9 


d+ + 


9 I - 9 I 


(X^ -|- QCi\ 


2fi'i5'2. 


cr_, 


( 12 ) 


where the effective frequencies and couplings are the same as in the previous case. At most, 
we may obtain a similar form to effective Hamiltonian Hid in the somewhat obvious regime 

gi 3> 2 g 2 and (72 3> 0 where we can approximate, 

2 

H2D ~ ^ lljdjSj -h A ^S|S2 S1S2 j -h gi ^Sj Si j d^,. ( 13 ) 

i=i 

Again, it is possible to interchange the helds via the transformation parameter Thus, 
Hamiltonian Hi will share the properties shown by H 2 in the particular regions gi S> 2g2 
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and g2 S> 2 gi. Furthermore, in the case of identical couplings, where H2 is invariant to 
rotations tj{ 6 ), 


H2D 


91=92 




2 

(cui + 0J2') 'y ^ a^jOjj + {0J2 
i=i 



\/2gi 


(32 ) cr+ + I ( 3 l 


(T_ 


( 14 ) 


This Hamiltonian is equivalent to two resonant helds, one of them interacting under the 
Jaynes-Cummings dynamics with the qubit and the other under anti-Jaynes-Cummings 
dynamics [ 18 ]. Note that under resonant helds, Ui = 0J2, it also conserves the quantity 
AC = —a\ai + a\a2 + -1-1/2 and, thus, can be solved. 


III. GROUND STATE CONFIGURATION. 

It is well known that the spin-iV /2 version of Hamiltonian H2 in the semi-classical limit. 
At 3> 1, allows four types of ground state na. These that can be characterized by four 
order parameters given by the qubit energy difference, ipz), the two mean photon numbers, 
(a^Oj), and the mean two mode photon number, (y) with y = (hi -1- 02). These 

conhgurations are dehned in Ref. [U] as: (i) a normal phase, where all order parameters 
are zero,((T^) = (ajoj) = (y) = 0, and corresponds to a separable ground state with zero 
excitation, (ii) two single-mode superradiant phases, where the qubit energy difference and 
one of the mean photon numbers are different from zero, (d^) 7^ 0 with (aihi) = (x) = 0 and 
(^2^2} 7^ 0 or (0302) = (x) = 0 and (aihi) 7^ 0, (hi) and a two-mode superradiant phase, 
where all the order parameters are nonzero, (a^) 7^ 0 , (djdj) 7^ 0 , and (x) 7^ 0 . For example, 
in Ref. [T 5 |, an exceptional solution is found under resonant helds, = LJ2 = cj, the two¬ 
mode superradiant phase appears just for equal couplings above a critical coupling, gi = 
92 > gc = \/wohi/ 2 . Here, we are going to discard the two-mode photon number operator as 
order parameter and use the more adequate two-mode SU{ 2 ) operator = (0102 -l- aia2)/2 
that describes hopping between modes as we already know that the systems show such a 
partial symmetry. 

In order to hnd the ground state conhguration for the single qubit models presented here, 
we set to study the exceptional case of resonant helds, Ui =002, that simplihes the problem 
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Figure 1 . (Color online) Four order parameters (a) mean energy difference, (d^), (b) first field 
mean photon number, (d|ai), (c) second field mean photon number, {a\a2) and (d) two-mode 
517 ( 2 ) mean hopping, (Jj,), for the ground state of Hamiltonian Hi on resonance, wq = ojj = u. 
The couplings are given in units of the field frequency, u). 


by eliminating the field-field coupling. 



(16) 

2 

= + a; ^ a]aj + g{a\ + hij 

j=i 

= a;d2a2 + Hr 

(16) 

(17) 


At this point, it is straightforward to recover the effective coupling transition for the quantum 
Rabi model in the weak coupling regime |T 9 l 1 ^ . 

9ic = ( 18 ) 

and conclude that the model presents a ground state configuration with no excitation, 
lfi')g|0)i|0)2; for coupling parameters in the range 0 < \Jgi+ gl < gic- This can be seen more 
clearly in Fig. [^a), where the mean qubit energy levels difference is shown for Hirf calcu¬ 
lated with standard numerical methods [2III22]. Figure [^b) and Fig. [^c) show the mean 
photon number in the hrst and second held, in that order, and Fig. [^d) shows the mean of 
the two-mode SU (2) hopping operator. It is also straightforward to borrow the deep-strong 
coupling, g ^ uq, result from the literature mm and realize that the ground state will be 
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(a) 




(^ 2 ) 


(b) (ni) 



1 



(c) (n 2 ) 


(d) (Jx) 



Figure 2 . (Color online) Same as Fig. for the ground state of H2 on resonance, wq = luj = io. 
The couplings are given in units of the field frequency, u. 


two-fold degenerate with ground state energy proportional to —g‘^/u, corresponding to the 
two parity separable ground states, (T|e) -l- \g)) |±/ 5 i)]^ |=t/92)2 with the fields in coherent 
states, 1 / 3 ) = with parameters / 3 j = gj/g. This ground state conhguration 

leads to mean values (d^) = 0 , (hi) = |/ 3 ip, (h2) = |/32p and {J^) = 2Re(/3i/32)- Note that, 
in the case of Hamiltonian Hi^ we can just describe two types of ground state configura¬ 
tions, the vacuum conhguration where all four order parameters are zero, and a non-vacuum 
conhguration where all four order parameters are diherent from zero. 

In order to hnd the critical coupling for the second model for resonant helds. 


H2RF = (a\ai + 0^02 j 

.t , 9 I- 9 I. ‘^9192 

g(i\ ~\~ -cii —- CI2 


9 


9 


d+ + 




9 


di -|- gtti 


‘^ 9 i 92 . 
-( 

9 


d_, ( 19 ) 


it is simpler to work with H2 in the zero and single excitation subspaces, then just half the 
critical coupling found [ 23 ] to obtain, 


5'2c — gic- (20) 

Thus, the vacuum ground conhguration will be in the parameter range 0 < \/gf + g^ < 92c 
and we recover the result in Ref. [I 5 ] for resonant helds and qubit. Again, this is simpler 
to see in the mean qubit population inversion. Fig. [^a). Figure shows the four order 
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parameters defined above for Hamiltonian H 2 . Note that we recover the four ground state 
conhgurations described in Ref. [15] if we just exchange their two-mode photon number, x, 
for the two-mode SU{2) hopping operator, J^. Thus, we are still able to see a change in 
the ground state conhguration for just the single qubit in the case of identical couplings, 
9i = 92 > 92 c, as expected from the semi-classical model analysis [Uj. We can write some 
of the different ground states for Hamiltonian H 2 - In the cases g 2 9c and gi > gc 
or gi gc and g 2 > gc the ground states will be given by (T|e) -l- \g)) lif?!)! | 0)2 or 

(=F|e) -|- \g)) |0);^ |± 5 ' 2)2 where the held states | ± gj)j are coherent states. In the regions 
gi 3> 2 g 2 with g 2 ^ 0 and g 2 2gi with (71 3> 0 the ground state conhguration will be given 
by that of Hamiltonian Hi in the deep strong coupling regime, (T|e) -l- \g)) |±/di)]^ |=t/^ 2)2 
with the coherent parameters as dehned beforehand. 

Note that, as expected, in these single qubit models the transition between ground state 
conhgurations is smooth and can be understood as a quantum precursor of the phase tran¬ 
sitions observed in the classical limit with an inhnitely large ensemble of qubits. 


IV. SPECTRA. 

Here, we will use continue our analysis of the exceptional case of resonant helds, ui = 
U 2 = OJ, to hnd a solution for the spectra of the models. In the case of Hamiltonian Hi, it is 
straightforward to see that the model is tractable in the exceptional case of resonant helds 
as in the displaced frame it can be written as Hirf- The proper basis for this Hamiltonian 
in the displaced frame is given by {|nb) 2 |±,thus, each an every eigenstate can be 
labeled by the displaced mean photon number of the second held, rih = 0,1,2,..., and 
parity component, (±, j) with j = 0,1,2,..., of the Rabi basis [3] EHT] [2ll |25|. Figure 
shows the hrst ten proper values of Hamiltonian Hi, note that each and everyone of them 
can be labeled at any given ehective coupling strength and the crossings in the spectra 
are always between diherent subspaces. The eigenvalues shown correspond to the following 
states, { 10 ) 2 !+, j)} with j = 0,1,2 in solid red lines, {| 0 ) 2 | —,j)} with j = 0,1,2 in solid 
blue lines, {| 1 ) 2 |-|-, j)} with j = 0,1 in dashed red lines, {| 1 ) 2 | —, j)} with j = 0,1 in dashed 
blue lines, {| 2 ) 2 |- 1 -, 0)} in dot-dashed red lines, and {| 2 ) 2 | —,0)} in dot-dashed blue lines. 
As expected for large values of the ehective coupling constant, g ^ u, the ground state 
will be twofold degenerate and the degeneracy of the rest will increase in twofold steps in 



Figure 3 . (Color online) First twelve members of the spectra of Hamiltonian Hi with branches 
labeled as follow: Positive parity branches, {^6)21+,^)} are shown in red, negative parity branches, 
{1^6)2! — ) j)}) are shown in blue. The branches corresponding to Ub = 0 , 1,2 are shown as solid, 
dashed and dot-dashed lines, in that order. 

the exceptional case of resonant fields; e.g., in this regime the two-fold degenerate ground 
state corresponds to {|0)2|±, 0)}, the four-fold first excited state, {|0)2|±, 1), |1)2|±, 0)}, the 
six-fold second excited state, {|0)2|±, 2), |1)2|±, 1), |2)2|d:, 0)}, and so on. 

As we shown before, Hamiltonian H2 has four regimes, {gi <C ca, (72 "C ca, gi ^ 2 g 2 , g2 S> 
2 gi}, where it can be approximated by Hamiltonian Hi. Thus, the spectra in these regimes 
can be labeled by the displaced second field photon number and Rabi basis as we have 
just shown above but there is a hfth regime gi = g2 where the spectra can be constructed. 
As mentioned in Sec. 2, the Hamiltonian H2 for the exceptional case of resonant helds 
and equal couplings conserves the excitation operator A/" = —a\ai + a\a2 + Ozl 2 -|- 1 . This 
partitions the whole Hilbert space in subspaces of inhnite dimension with the same mean 
value {M) = 0 , ± 1 , ± 2 ,.... Note that we have the three labels we need to uniquely identify 
each eigenstate, {\±,nd, j)}, the total parity of the state, ( tt ) = ±, the displaced excitation 
operator dehning the subspace, Ud = (A/"), and the ordering in the subspace, j. Figure 
1^ shows six members of the spectra following this convention, the ground state will always 
be |-|-,0,0) and will be completely degenerate at large effective coupling values gi ^ oj 
due to the fact that each subspace has an inhnite dimension. In Fig. |^we also show the 
eigevalues corresponding to |-1-, 0, j) with j = 0,1 and |-1-, 2, 0) in solid and dashed red lines, 
respectively, as well as those related to |—,l,j) with j = 0,1 and | —,—1,0) in solid and 
dashed blue lines, in that order. 
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Figure 4 . (Color online) Six members of the spectra of Hamiltonian Hi with branches labeled as 
follow: The eigenvalues corresponding to |+, 0 ,j) with j = 0,1 and |+, 2 , 0 ) are shown in solid and 
dashed red, respectively, those related to |—, 1, j) with j = 0,1 and | —, —1, 0) in solid and dashed 
blue, in that order. 

V. SINGLE EXCITATION DYNAMICS 

Dynamics of the quantum Rabi model are well studied in the most relevant regimes; weak 
[261 EH, g ^ w, ultra-strong | 28 ], g > O.lca, and deep-strong coupling [ 29 ] , g > uj, regimes 
[ 29 ] . Here, as a practical example, we consider the exceptional case of resonant helds in order 
to entangle the hrst and second cavities with a single excitation near the deep-strong coupling 
regime. In the weak-coupling regime, g w, the time evolution operator for the resonant 
quantum-Rabi model is well known. In the single excitation subspace, {|e)q| 0 )i, |5f)q|l)i}, it 
yields the following time evolution for an initial state |t/’(0)) = |e)q|0)i|0)2, 

^ cos5ft|e)g|0)i|0)2 -ismgt\g)q [cos^|l)i|0)2 -b sin^|0)i|l)2], g <^uj. (21) 

The time evolved state oscillates between the original state and an entangled state of the 
two helds. Thus, for a single excitation in the weak coupling regime, the hrst model can 
act as a SU( 2 ) beam splitter at times t = n 7 r/{ 2 g). There, the qubit will be at the ground 
state and the cavities will be sharing a photon in a ratio | cos.^p/| sin,^p = gf/g^ given 
by the ratio between the two qubit-held couplings as tan^ = (72/5'! • Now, in the weak 
coupling regime, this process is slow, in order to obtain a fast beam splitter we need to 
go for stronger couplings. Figure [^ shows the time evolution of the mean qubit population 
inversion and mean photon numbers of the helds for a 50/50 beam splitter realization in the 
ultra-strong regime, gi = g2 = 0 . 15 a; that gives g"^ = 0.212 o;^. The numerical evolution is 
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Figure 5 . (Color online) Time evolution of the (a) mean energy difference, {d'z), (b) first field 
mean photon number, (a|ai), (c) second field mean photon number, (0202) and (d) hdelity, T" = 
|('0(t)|e“*^i*|V'(O))p, for an initial state |V'(0)) = |e)g|0)i|0)2 under dynamics given by Hamiltonian 
Hi on resonance, uq = loj = lo and equal qubit-field couplings in the ultra-strong coupling regime, 
9 i = 92 = 0 . 15 a; that realizes a 50/50 beam splitter under weak coupling. 

compared to the result obtained in the weak coupling regime and it is possible to see that 
they are in close agreement during the first oscillation. If we stay with this effective coupling 
regime, = 0.212 (o;^, we can keep this level of agreement for different splitting parameters. 
Fig. @ presents the case g\ = that in the weak coupling regime realizes a beam 

splitter. In both cases we can see that the qubit does not reach complete transfer to the 
ground state, still a conditional measurement of the qubit in the ground state delivers a 
state close enough to the ideal split single photon state; Figures [^d) and |^d) show the 
hdelity between the weak coupling evolution and the exact numerical ultra-strong coupling 
evolution, IF = |(-?/i(f)|e“*^i*|V’( 0 ))p. 

We can do the equivalent with Hamiltonian H2 [301 E], but in this case the one- 
excitation subspace for weak couplings, gj ca, will be dehned by the tripartite basis 
{|5f)q|l)i|0)2, |e)q|0)i|0)2, |5')g10) 111)2}, and obtain the evolution for an initial state |t/’(0)) = 
|e)q|0)i|0)2, up to an overall phase constant, 

|Vi2(f)) ~ cos5(t|e)q|0)i|0)2 -ismgt\g)g [cos^|l)i|0)2 - sin^|0)i|l)2], g <^u. (22) 

Again, in the weak coupling limit, the time evolution under H2 dynamics delivers a single¬ 
photon beam splitter state of the two helds at the time t = mi/{ 2 g). Figure and Fig. 
show the numerical exact evolution of the initial state in the ultra-strong coupling regime 
compared to the approximate result obtained using the weak coupling evolution for a 50/50 
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Figure 6. (Color online) Time evolution of the (a) mean energy difference, (b) first field 

mean photon number, (a|ai), (c) second field mean photon number, (0202) and (d) hdelity, T = 
|('0(t)|e“*^i*|V'(O))p, for an initial state |V'(0)) = |e)g|0)i|0)2 under dynamics given by Hamiltonian 
H\ on resonance, uq = loj = lo and equal qubit-field couplings in the ultra-strong coupling regime, 
gf = 3^1 with identical effective coupling to that in Fig. = 0 . 212 a;^, that realizes a 75/25 

beam splitter under weak coupling. 



Figure 7 . (Color online) Same as Fig. [funder dynamics given by Hamiltonian H2. 


and 75/25 realization with the same parameters as those in Fig. |^and Fig. in that order. 
Note that the hdelity in these cases is better than those nnder Hamiltonian Hi dynamics 
for short evolntion times bnt it seems to degrade faster for longer evolntion times. 


VI. CONCLUSIONS 

We have stndied exceptional solntions for two models describing a single two-level system 
conpled to two boson held modes. The hrst is a hnite dissipative two-state system with only 


12 

























Figure 8. (Color online) Same as Fig. [funder dynamics given by Hamiltonian H2. 

two fields in the absence of tunneling. The second is equivalent to two-orthogonal helds 
coupling to the two corresponding orthogonal dipoles of a two-level system. Both models 
may be feasible of experimental realization in cavity- or circuit-QED da eg. 

The models conserve parity and show a partial SU( 2 ) symmetry involving the two boson 
modes, thus, we explored regimes where they may be related to well known models with 
similar structure, like the quantum Rabi model. We focused on the exceptional case of 
resonant helds where the models are analytically tractable. Although only one of the models 
can be transformed to a form including the quantum Rabi model, we found that the ground 
state conhgurations of both models present the same critical coupling than the quantum 
Rabi model. Around this critical coupling, the ground state goes from the so-called normal 
conhguration with no excitation, the qubit in the ground state and the helds in the quantum 
vacuum state, to a ground state with excitations, the qubit in a superposition of ground 
and excited state while the helds are not in the vacuum anymore, for the hrst model. The 
second model shows a more complex ground state conhguration landscape where we hnd the 
normal conhguration just mentioned before, two single-mode conhgurations, where just one 
of the helds and the qubit are excited, and a dual-mode conhguration, where both helds and 
the qubit are excited. For the hrst model and some regions of the second model, we showed 
that the held components of the ground state are given by coherent states in the deep-strong 
coupling regime. Following the integrability criteria for the quantum Rabi model established 
by Braak [ 3 ], our results point to the integrability of these exceptional solutions in these 
regions. 

We have also shown that these models for resonant helds and qubit frequencies can be 
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used as two-mode SU{2) beam splitters in the case of single excitation even in the ultra¬ 
strong coupling regime; at specihc times, an excited two-level system may give its energy 
quanta to the two commuting boson helds, this excitation will be shared in an entangled 
state by the two fields in a manner dictated by their coupling ratio. 
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